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Abstract. We evaluate explicitly, in terms of the Cauchy data, the constant pre-factor in the 
large x asymptotics of the Painleve III tan-function. Our result proves the conjectural formula 
for this pre-factor obtained recently by O. Lisovyy, Y. Tykhyy, and the hrst co-author with 
the help of the recently discovered connection of the Painleve tan-functions with the Virasoro 
conformal blocks. Our approach does not use this connection, and it is based on the Riemann- 
Hilbert method. 


1 Introduction 


We consider the particular case of the third Painleve equation, which is a radial-symmetric 
reduction of the elliptic sine-Gordon equation 

Ux 

Uxx H-1- sm u = 0 

X 

Starting from the pioneering works n, na on the Ising model, this equation has been playing 
an increasingly important role, as a “nonlinear Bessel function”, in a growing number of physical 
applications (see e.g., |1] and references therein). Apparently, the hrst appearance of equation 
([^ in the physical applications should be credited to work of J. M. Myers [TU] . 

Equation ([^ can be written as a (non-autonomous) Hamiltonian system. 


( 1 ) 


du dpi dv dPL 

dx dv ’ dx du' 
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on the phase space = {('U,n)} equipped with the canonical symplectic structure, 
fl = dv A du. 


( 2 ) 


The Hamiltonian "H is given by the formula 


H =- X cosu. 

2x 


We are concerned with the global asymptotic analysis of the r-function corresponding to the 
Painleve III equation Q which is dehned according to the equation (see [H], [TU]h 


d In r 1 ^ , 

—j — = “7^ 

dx 4 


(3) 


In fact, it is this tan-function (evaluated for a special family of solutions of equation ([^) that 
played a key role in the above mentioned Barouch-McCoy-Tracy-Wu theory, and, since then, it 
has appeared in many problems of statistical mechanics and quantum held theory. Let us now 
remind some of the basic known facts about the asymptotics of the solutions of equation ([^. 
We refer the reader to monograph [1] for more details and for the history of the question. 

Equation Q possesses a two-parameter family of solutions characterized by the following 
behavior at a: = 0, 


u{x) = a\nx + (3 + O ^ a; —)■ 0, 


(4) 


where the complex numbers a G C, IQ'al < 2 and /3 G C can be taken as parameters - the 
Cauchy data, of the solution u{x) = u{x\a,(3). The behavior of the solution u{x\a,f3) as 
X -A -l-cxD, is known. For an open set in the space of parameters a, f3, which we will describe 
later, the large x behavior of u{x\a, (3) is oscillatory, and it is given by the formulae. 


u{x) = [ - ]] + b_e 


-“a:-*‘^-i/2 1 + O - + 


-1-0 (mod27r), a; —)■ oo. 


(5) 


where 

i/ = —^6+6-, |Q^z/| < 1/2. (6) 

The asymptotic parameters at inhnity - the complex amplitudes b±, can be, in fact, expressed in 
terms of the Cauchy data a, f3, and the condition IQ^z/l <1/2 imposes additional restriction on 
them. The corresponding connection formulae were obtained in 1985 by V. Yu. Novokshenov 
[18] (see also: [101 and m), and they are given by the equations. 


= 


sin 2717] 
sin 27ra ’ 


where 


6, = _e¥Tf 2‘±-ir(l 




sin 2777] ’ 


1 i 

‘"■^4 + 8“’ 


f /3 -|- a In 8^ — In —- tT 

dvrV^ ) 277 r(4 + f)’ 


(7) 

( 8 ) 
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and r(z) is Euler’s Gamma-function. The open set in the space of the Cauchy data a, (3 where 
the both asymptotics, Q and ([^ are valid is described by the inequalities (see also Remark 
below), 


0 < 3fJcr < - 


IQ'al < 2, sin 27iri ^ 0, 


sin 27r?7 
— 

sm27r(T 



< 2- 


(9) 


where a and r] are understood as functions of a and [3 dehned in ([^. We notice that this set 
contains all sufficiently small pairs {a,/3), all real pairs (a,/?) such that the corresponding rj 
satishes the inequality 0 < ^ (mod (1)) and all pure imaginary pairs such that |a| < 2. In 

fact, it is convenient to take a and rj as the independent parameters and think about a and f3 
as their functions, i.e., 

a = i{2 — Sa), /? = —tt- f dvrr/— z(2 — 8cr) ln8 — 2i In —(10) 

r(2cT) 

where a, rj are the complex numbers satisfying ([^. The expressions of the asymptotic param¬ 
eters at X = cxD in terms of a and r] have already been presented in ([^. 

The derivation of formulae ([^ is based on the Isomonodromy-Riemann-Hilbert Method. We 
again refer the reader to monograph |1] for more details and for general references concerning the 
connection problem for Painleve equations. In the framework of the Riemann-Hilbert method, 
the parameters a and r] have an independent important meaning as the monodromy data of 
the auxiliary linear system associated with the third Painleve equation. This meaning of the 
parameters a and rj plays important role in the considerations of this paper, and it will be 
explained in detail in the next section. 

Equations (|^ and ([^ in turn imply the following behavior at zero and at inhnity of the 
corresponding tau-function (see also 0 ), 

r(x) = Gox“^(lo(l)), X-)■ 0, (11) 


and 

t { x ) = + 0 { 1 )), X —)■ CXD. 


( 12 ) 


In fact, one can write a complete asymptotic series for the tau-function at both critical points 
whose coefficients are explicit functions of the Cauchy data a, (3 or, equivalently, of the mon¬ 
odromy data (j, r]. The issue which we are concerned with is the evaluation of the ratio 

Goo/Go (13) 


in terms of the initial data a, (3. This can not be done just by using the asymptotic equations 
(|^ - ([^ and the connection formulae ([^-(|^. Indeed, we are dealing here with the “constant of 
integration” problem. For the special one-parameter family of solutions of equation ([^ related 
to the Ising model, this problem was solved by C. Tracy [20] in 1991. This special family is 
obtained by putting 


V 


0 and cr G M, 


0 < (T < 


1 

2‘ 


(14) 
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in (10). Zero value of r] is excluded from set ([^ which means that the behavior of this special 
family at inhnity is very different from the oscillatory one given in ([^. In fact, all the solutions 
from this family exponentially approach vr (mod27r), 


u{x) — vr ~ iK\j^x ^^ 


a; —)■ oo, K = —2cos27rcr. 


In his calculations, Tracy made use of the existence of an additional Fredholm determinant 


representation of the tan-function evaluated on the family (14). We are interested in a generic. 


two-parameter case where there is no such representation. A conjectural answer to the problem 
has been produced in [T2] with the help of the recently discovered by O. Gamayun, N. lorgov, 
and O. Lisovyy connection of the Painleve tan-functions with the Virasoro conformal blocks 
0.0- In this paper we prove the conjecture of na. Our main result is the following theorem. 

Theorem 1. Let a and rj be the “monodromy” parameters of the Painleve III function u{x) 
satisfying the inequalities 


Then the ratio (13) is given by the formula, 


Co 


226 - 






r(l — 2<j ) /G(l + w)G(l + 2(7)G(1 — 2(7 )G(1 + (7 + // + — a — ti)\ 

r( 2 a) V G{l + a + ri+ ^)GC-^ - (t - p) 7 ’ 

where v is defined in 0 and G{z) is the Barnes G - function. 

It should be noticed that in [12] a slightly different dehnition of the tan-function is used. 


The exact relation of the constant (15) and the one conjectured in [12] is discussed in the last 
section of the paper. 

Our proof of Theorem [T] is not based on the conformal block connection. We use the 
Riemann-Hilbert representation of the third Painleve transcendent and the Malgrange-Bertola 
extension of the Jimbo-Miwa-Ueno dehnition of the tau-function. 


In the course of our proof, we also conhrm one of the key observations of ua that ratio ( [15| 
determines the generating function of the canonical transformation of the canonical variables 
determined by the initial data {a, fi) to the canonical variables determined by the asymptotic 
data (6+, &_) (see the end of Section 5 for more detail). In fact, this Hamiltonian interpretation 
of the pre-factors in the asymptotics of the Painleve tan-functions was hrst suggested in the 
work [7] of N. lorgov, O. Lisovyy and Yu. Tykhyy. 


The evaluation of ratio (13), which we have made rigorous in this paper, is only one of 


a series of highly nontrivial predictions and already established facts which came from the 
remarkable discovery of Gamayun, lorgov, and Lisovyy. These other predictions and results, 
including the key ingredient of the approach of [3] and [3] , which is the explicit conformal block 
series representations for the Painleve tau-functions, do not yet have their understanding in 
the framework of the Riemann-Hilbert method. 

We shall start the proof of Theorem with the reminding of the Isomonodromy-Riemann- 
Hilbert formalism for the Painleve equation ([^ (for more detail see, e.g., [1]). 
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2 The Riemann-Hilbert Representation of the Solutions 
of the Sine-Gordon/Painleve III Equation 

The Riemann-Hilbert problem associated with equation Q is posed on the oriented contour 
r depicted in Figure]^ and it consists in the hnding of a 2 x 2 matrix-valued function T(A) 
which satishes the following properties. 



Figure 1: Contour F 


• The function \1/(A) is analytic on C \ {F}, it has continuous ± - limits on the contour 
F, and these limits satisfy jump condition \1/+(A) = \1/_(A)S'(A). Here “ -f ” denotes the 
boundary values from the left side of the contour and “ — ” denotes the boundary values 
from the right side of the contour. Jump matrix S{X) is piecewise constant, its different 
components are indicated in Figure [TJ and they are given by the equations, 


^(oo) _ ^( 0 ) 


1 OA 

p + q 1 / ’ 


^(oo) _ ^(0) 


1 p + q \ 

0 1 ’ 


(16) 


E = 


a/1 +pq V 


P 

1 


p,qeC, 1+pqj^O, ai = 


(17) 


• The function T(A) satishes the following conditions at zero and inhnity 
T(A) = Po{I + MfU + A ^ 0, 


^(A) 



A 


+ 0 



-C73 


A —y oo, 


( 18 ) 
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where Pq, M^\ here are some constant in A matrices and 


c^s = 


1 0 
0 -1 


The Riemann-Hilbert problem is nniqnely and meromorphically in x solvable for all 
p,q E C, 1 + pq ^ 0 [IT] and the corresponding solntion u{x) = u{x;p,q) of the third 
Painleve eqnation ([^ is given by the formnla, 

u{x) = 2 arccos(Po)ii- 

In fact, the following eqnation takes place, 


P 


0 


cos(|) 

-isin(|) 


-isin(|) 


cos(|^ 


= e 




(19) 


Remark 1. The above Riemann-Hilbert setting corresponds to the generic solutions of &■ 
There is one parameter family of a separatrix solution which is characterized by the following 
Riemann-Hilbert data 


5*1°°) = = 


1 0 
K 1 




1 K 
0 1 


E = ±i 


0 1 
1 0 


K e 


This is the family which includes the McCoy-Tracy-Wu solution (If) and which is not considered 
in this paper. As it has already been mentioned, the constant problem for this family was solved 
in 


The parameters p,q G C in (16), (17) are connected with the parameters of asymptotic of 
u{x) via 

.sin 27r(cT + p) .sin 27r(cT — p) 


p := 


q := t- 


sin 2np sin 27rp 

Conditions ITJal < 2 and l^yz/l <1/2 can be rewritten in terms of p and q as 

pq ^ {—ioo,—2i] U [2i,-\-ioo), and pq ^ {—oo,—1], 

respectively. 


( 20 ) 


Remark 2. Conditions (20) are the conditions which appear during the asymptotic analysis of 
the Riemann-Hilbert problem. The asymptotic parameter v is related to p, q according to the 
equation, 


10 = ln(l +pg). 


Also, 


l+pq = 


sin^ 2Tia 
sin^ 27r?7 


We restrict ourselves in ^ to the inequality |arg \ < vr/2, instead of the inequality 
arg by a technical reason. This means that we actually analyze one of the 
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components of the full set of the Cauchy data corresponding to the generic asymptotic behavior 
0 and 0. The another component is defined by the condition, 


sin 2TTri 
arg—— 

Sm/TTCT 



<7r/2, 


which in turn implies the following change in formulae 0, 

^ sin 2irr; 
sin 27r(j 

The analysis presented in this paper can be easily extended on this component of initial data 
as well. Actually, the t- function does not change if we add 27ii to the function u{x). But rj is 
shifted by So such change of variable allows us to go from one component to another. 

Function tl/(A) satisfies system of linear ordinary differential equations 


^ = /1(A)>I>(A), (21) 

d'^ 

- = mmx,. 


A(X) - ^ Po(*0-3)-Po ^ 


U{\) = 


16 4A 
iXxa^ iUx<Ji 


A2 


( 22 ) 


Equation Q is the compatibility condition for this system and it describes isomonodromic 
deformations of the system (21). From this point of view a and rj play role of the monodromy 
data. 

We complete this overview of the Riemann-Hilbert formalism for equation ([^ by presenting 
the general alternative definition of the Jimbo-Miwa-Ueno tan-function in terms of the solution 
4/(A) of the Riemann-Hilbert problem. 

Dehne ^ 

ij(oo)(A) ^(A)e^-3^ |A|>P. 

Then, according to [9] the equation, 

WJAW = - res Tr((#<“>(A))-‘(#<“>(A))'(-^<73))<li, (23) 

A=oo \ O / 

defines the differential form whose antiderivative is the logarithm of tan-function. Actually, 


from (18) we have that 


^F(“)(A) = / + 


Ml 


(oo) ^ 

A +°(a5 


(24) 


Substituting (18) to the equation (21), one can express in terms of u{x) (see |4].[T7]) 


= 


2ixu.r 


x^ 


-(72 — I COSU 


u: 


O's, <72 — 


0 -i 
i 0 


(25) 
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Substituting, in turn, (24)-(25) into (23) yields the relation, ojjmu = dlnr(a;), where dlnr is 
dehned in In other words, the tau-function can be alternatively dehned as 


r = TjMu{x,p,q) = 

So dehned the tau-function is unique up to multiplication by a constant depending on p and q. 
The key fact for us is that, following [2], it is possible to extend Jimbo-Miwa-Ueno differential 


form (23) on vector helds in p and q in such a way, that it will remains a closed form. Such 


extension will allow us to dehne tau-function already up to a constant, which does not depend 
on p and q. 


3 Malgrange-Bertola Differential Form 

In this section we basically repeat the calculations and the results of Section 5.1 of paper |2] 
adjusting them to our special case. 

Put F(A) = Denote G(A) the jump matrix for D(A). Following |2],[I1], we 

dehne the Malgrange-Bertola differential form by the equation 

uJMB[d] = j Tr(y-ir (26) 

r 

Here d denotes the vector held in the space of parameters a;,p, g, and the prime denotes deriva¬ 
tive with respect to A. This diherential form was introduced originally by B. Malgrange in HI 
for the case when the contour T is a circle. M. Bertola in [2] has extended the Malgrange’s 
dehnition to an arbitrary Riemann-Hilbert setting. 

Let us establish the connection of this form with Jimbo-Miwa-Ueno form. For the case of 
general Riemann-Hilbert problem, the analog of this Lemma was proven in j2]. 

Lemma 1. The Malgrange-Bertola differential form, evaluated on the vector fields in parameter 
X, is equal to the Jimbo-Miwa-Ueno form up to a term, depending only on G(A). 




8 ^^J2TTi 


X 

4' 


(Xmb — x}jMU — / Tr 

Proof. First, we have 

G(A) = 

where S'(A) is the jump matrix for J/(A). Hence, 

dxG{X)G ^(A) = J- ——G{X)a3G ^(A). 

o o 

Also, we have that 

y+(A) = y_(A)G(A), 


(27) 


(28) 


(29) 


and, by A - diherentiation, 

U[(A) = u:(A)G'(A) + U_(A)G"(A). 


(30) 









Substituting (28),(29),(30) in (26) we get 


^MB[dx] — 


j iv((y->y; 


Y-^Y')—a, 
- 8 




27ri 


Let us introduce the notation for the parts of the contour L as it is indicated in Figure 
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We have 


r5ur6ur4 


Ta 


Tsurg 


We also have 
Hence, 


r(A) = |A|>i?. 


lAYrV'^a,']^ 


V + + 8 V2m J \ 

riuFs riur2ur3ur4 




8 7 2Tii 


A=oo \ 8 8 

,ixX 


= re^Tr((il“)(A))-'(4''“>(A))'^<T3) - j, 


and (27) follows. 


□ 


This lemma means, that Malgrange-Bertola form is indeed a good candidate for extension 
of Jimbo-Miwa-Ueno form. However, there is an additional term, depending only on G(A). One 
can, following again |2], cancel it considering the modified Malgrange-Bertola form u = UMB + d, 
where 

e\B] = 1 1 1V(G'G-‘(SG)G-')|^. 

f 

In the notations of [2] a; is the form O from Dehnition 2.2 of [2]. 

We have u:[dfi\ = u:jMu[dx] - f • Indeed, 

Tr(G'G-'(a.G)G-i) = ^ 2Tr(/ - G-^a,Ga,), 

and the additional term depending only on G(A) cancels. 

In the next section we will express the form u in terms of the coefficients of the asymptotic 
expansions of H(A) at A = 0 and at A = oo. We call this expression a “localization” of 


the original integral formula (26) for the Malgrange-Bertola form. This localized version will 


simplify dramatically the further analysis of the form uj. 


4 Localization 

Let us introduce the function 


0(A) = 5F(A)y(A)-\ 
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where d means the differentiation with respect to one of the three parameters, x, p, or q. The 


5-version of equation (30) reads, 

ar+(A) = aF_(A)G'(A) + F_(A)aG(A). 


Expressing G(A) from ( [2^ as G(A) = Y_ ^(A)y'+(A) we rewrite the last equations as 

ay+(A)y-i(A) = aF_(A)yri(A) + y_(A)aG(A)y-i(A), 


or as 


aF+(A)F+-^(A) - aF_(A)yr^(A) = F_(A)aG'(A)F+-^(A). 

The Sokhotski-Plemelj formula would then imply (cf. Lemma 2.1 of j2]) that 
Y_{y)dG{y)Y^\y) dy 


0(A) = 


y-x 


27ri 


Substituting E(A) = \k(A)e*^*'i6 (21), we have 


y'(A) = Ai(A)y(A) + (^ - ^)y(A)ff3, 


and 


.MslO] = / T.(y_-r OG)G-)^ = / (^ - ^)T.(.3(aG)G-)^ 


_i dX 


(31) 


(32) 


(33) 


+ f Tr(AY.(dG)Y-')^^. 


We introduce notation 

y(A) = Po(-r+mi A + 0(A2)), a^o, 


(34) 


(36) 




y(A)- (/ + ^ + o(^)), A^oq. 

Substituting these expressions for Y (A) to the definition of 0(A) we get 
0(A) = idPo)Po-^ + Poid mi)Po-'A + 0(A2), A ^ 0, 


0(A) = 


9m 


(oo) 




Comparing these formulae with (32) we arrive at the relations 


Y-(y)dG(y)Y^ Hv)^. = 


(36) 


(37) 
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Y_{y)dG{y)Y^ {y) dy 


y 


2Tii 


: = (SP(,)P„-‘, 


(38) 


Y4y)dG[y)Yi'[y) dy 


y 


2Tii 


; = Po(9 mi)P„ 


(39) 


Let us look now at the last integral in equation (34). Putting in it formula (22) for ^(A), we 
will see that this integral can be re-written as 

r r 


ZXUn 


Tr di 






““'.T. I .. / + T. f P,(..3)n- ^ 


A 27ri 


A^ 27rz 


The last equation, with the help of (|37|) - (|39|), is transformed into the localized formula, 

^-TT{a,{dPo)P,-^)+tTT(a3dm,y (40) 


Tr{AY4dG)Y-') ^ = ?TTv(a3Sm<“') - 


Substituting the derivative of G with respect to A in the formula for 0, we have that 


Together with (40) this gives us the following formula for u 

.2 




ix^ 

-f r 

16 


f (^a3dm[°°^^ - ^^^Tr (ai(dPo)Pg -h iTr (^asd ruij . 


One can check directly that 


Tr(^cT3((aG)G-i + G-\dG))^ = 0. 
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Therefore, we finally have that 


= ^Tr j - ^^^Tr (ai(aPo)^o + ^Tr {a2,d mi 


(41) 


Substituting (35), (36) to the equation (33), one can express the coefficients of asymptotics of 


F(A) at A = 0 and at A = oo in terms of u. In particular, one gets (cf. [T7]i. 


(oo) 

m; = 


2ixu. 




u. 


(72 — i cosM — 1 —<^ 3 ) 


o IXUx ■ ( 

mi= -^(T2 - I ( Y^{cosu - 1 ) - 


x'^ul 


era- 


Inserting these equations together with formula (19) in (41) and using also the fact, that u{x) 


satishes ([^ we transform equation (41) into the hnal expression for the form ui in terms of u 
and its derivatives with respect to x, p and q. 

Proposition 1. The modified Malgrange-Bertola differential form u admits the following rep¬ 
resentation 


u = 


xui X. 7 { xUxUp . , 

-h — (cosM — Ij I dx — I -—UpSmu + -—UxUpx H - ) dp 

8 4 / V 4 4 4 


2 2 

X " 'y'g I ■ 7 / 

—Mg sm u + -^UxUqx H-^ ) dq. 


(42) 


We notice that from (42) we have again the statement of Lemma that is that ijj[dx] = 
cla;lnr — We want to mention again that this part of the localization formulae has already 
been obtained in |2]. We also want to emphasize the important role which is played by the 
A - equation (21) in the derivation of the p, g - part of equation (42). It is the use of this 


equation that allowed us to present the original Malgrange-Bertola integral (26), hrst in the 


form (34), and then in the localized form (40). In fact, similar technique has already been used 


in the study of Toeplitz determinants with the Fisher-Hart wig singularities in paper |3] 
Appendix 6 and Lemma 6.2 of that paper. 


see 


Remark 3. As it was pointed out to the authors by M. Bertola, equation (31) can he used 
in the derivations of this section one more time and help to make a significant short cut from 


equation (34) to the localized form (fO). Indeed, Bertola’s suggestion is to use relation (31) for 
the product Y'_(A)G(A)F,r^(A) in the last integral of (34) directly and rewrite this integral as 


TV(/1(A)F_(A)5G(A)F-‘(A))P = Jti 


a1(A)(sF+(A)F-'(A) -SF_(A)F_-'(A))) (f 


resTr(^A{\)dY{X)Y-\X)y (43) 

poles of A{X)dX 


This is a quite general construction which allows one to localize the Malgrange-Bertola form for 
an arbitrary isomonodromic Riemann-Hilbert problem. In our case, one has to evaluate the, 
properly understood, residues at the points X = 0,oo. The result will be equation (fO). 
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5 Proof of Theorem [U 

Let us compute du. First we have, 


d 


^ / IN \ > 

— -— (cosM — 1) ) ax 

8 4 


XUxUr)x XUnSmU\ , , f XUxUax XUa^UlU . , 

^ - 1 dp^dx+ ^ ^ - dq^dx. 


4 4 y " V 4 

Then, using the fact that u{x) satishes equation ([^, we get that 

d 


xUxUjjx xu„smu . , 

px ^ - ]clx A dp 


2 2 

ry>^ 'T*'? / ■? / 

vL vL X 1 7 

—Up sin M + —u^Up^ H-^ ) dp 

2 2 2 2 

tX^ tX^ tX^ tX^ 'X 00 

-j- I —^VjpqSlYiU ~\~ ^ VjpUq COS ZA “1“ ^ ^px^qx 4“ ^ ^x^xpq ~t“ 4“ ~^^p^qx ) ^ dp. 


and 


d 


2 2 

/-y.^ ry>^ I 01 

Jj ^ Jj dj LLx (Xq . 

—MgSmu + —u^Ug^ H- - — ) dq 


XUxUax XUaSmU . , 

^ + —4_- ) (io: A cfg 


X^ . X^ X^ X^ X x\ 

T I ^ sin U T ^ UpUq cos U T ^ UpxUqx ^ UxUxpq T ~^^xUpq T ~^UqUpx | dp A 




Adding up the last three equations we obtain that 

VpUq Vqtlp 


du = 


dq A dp, 


where n = xUx- From equation ([^ it follows that 


d 

(z^pZ/g- VqXlpj 0, 


and hence we can observe that 


duj = lim duj = 

a:^0 


(y.pPq (y.q[5p 


dp A da 

, dq A dp = -:-. 

4 1 ^ ^ 


Therefore, if we dehne 


X adp 

w = u + —ax H-—, 

4 4 


then the form w will be a closed form on the full set of parameters, (x, p, q) and such that 
w[dP\ = wjMu[dP- This means we can put, 

r = 


and this equation would dehne the tan-function up to a constant, which does not depend on p 
and q. 

Remark 4. It is worth noticing, that from our analysis it follows that, in the case of the 
Painleve III equation 0 , the external differential of the (modified) Malgrange-Bertola form 
oj is proportional to the canonical symplectic form ^ for the Hamiltonian dynamics of the 
Painleve equation; indeed, we have that i.e.. 


du = —f2. 
4 
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The next step is to evaluate the small and the large x asymptotics of the form w. To 
this end we shall use asymptotics (|^,([^ for u{x) and make the following temporary technical 
assumptions, 

< 1. (44) 

6 

In our calculations we will need more terms of the large x asymptotics at inhnity which are 
given in 

u{x) = 

+ —(6z/2 + Aw - - Aw - 

8 8 


^ ’,3 ^Zix ^ j,3 „—3ix^—3iu—^ i ^—^+5\Qu\ 


■—bie^^^x 


48 ■ 


+ 0{x~ 


)■ 


(45) 


Substituting this asymptotics at the right hand side of equation (42), we shall arrive, after 


rather tedious though straightforward calculations, at the following asymptotic representation 
of the form a; as x —>■ cx). 


uj = d{2ux + z/^lnx + z/^) - -{b+db_ - b_db+) + ( —x 


+0(x-2+6l^^l)dx + 0(x-^+®l^^l))dp + 0(x-^+®l^^l))dg, X ^ CX), 


(46) 


The derivation of the small x asymptotics of the form u is based just on the estimate (|^, i.e., 
no need for its extension, and it is much easy to obtain. 


, , adfd 

u = d [ -Inx-— 


+ 0(x 


l-|S(a)| 


))dx 


+0(x^ lnx))(ip + 0(x^ Inx))^^, x —)■ 0. 


(47) 


As it has already been indicated, the derivations of formulae (46) and (47) are straightforward. 
However, because of the importance of these formulae for our further analysis, we present the 
details of their derivations in the Appendix. 

In view of the assumptions (44), estimates (46) and (47) yield the following asymptotic 
representation for the form w, 


o? o?' 

w = —d ( — Inx + — ) + o(l), X —)■ 0, 


(48) 


and 


w = d (2z/x + Inx + z/^) — ^{b+db_ — b_db+) + ^dx + + o(l). 


X —)■ CX). 


(49) 
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On the other hand, from ([IT| and (12) we have that 


a 


w = —d I — Inx ) + dlnOo + o(l), x —)• 0, 


and 


X 


w = d \ 2ux + v^hix + — + din Coo + o(l), x ^ oo. 


The comparison of (48) - (49) and (50) - (51) implies that 

>, 2 ' 


din Co = —d I — 


a 


and 


din (Too = di/^ — -(6+d6_ — b-db+) + 


ad/9 


(60) 


(51) 


(52) 


The last two equations mean that 


(Too , / o a^ . \ adjd i 


(where we have also taken into account (|^), or that 

In = + ~ 2ibj^dbJ) + c, 

Co 8 4 T 

where c is the numerical constant, independent on p and q. 
Following [12], we introduce notation 

g-4«p ^ sin27r(a + r7) 
sin 2ttp 


(53) 


(54) 


Using this and the connection formulae ([^, (10), we can re-write the differential form 
\{(y.d(3 — 2ib+db_) as the differential form in variables p, a and i/. 


-(ad/9 — 2ibj^dbJ) = —Sni{adp + ivdp) -|- 2nidp — (12 — 48(t) ln2d(T 


-t-(fvr -|- 41n2)i/dz/ -|- (1 — 4(T)dln ^ ^ -|- 2iz/dlnr(l -|- iu). 


T{2a) 


Therefore, we can re-write (54) as 


C 

In -|- 4cr — 8 ct^ — iv ^ 27rip — 12a In 2 -|- 24 ct^ In 2 -|-h 2z/^ In 2 

(To 2 


-8ni J{adp + iudp) + j (1 —4a)dln ^ ^ 4- j 2ivdhiT{l + iv) + c. 


(55) 
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It remains to evaluate the integrals in (55). For the integrals involving the F-functions one gets, 


(1 — 4cT)(iln 


F(l-2a) , F(l-2a) 


r(2<T) 


= In 


r(2<T) 


4a + 8a^ + 2ln(G(l-2a)G{l + 2a)]+c, (56) 


J 2ivd\\iT{l + iv) = iv — — zz/ln(27r) + 21nG'(l + iu) + c, 

where G{z) is the Barnes G-function and we have used the classical formula, 
[ lnF(a;)(ia; = ^ ln(27r) + 2 ;lnF(; 2 ) — lnG(l + z) -\r c. 

f Li Zi 


(57) 


The most challenging, i.e., the hrst integral in (55) has already been evaluated in n. Here is 
the result. 


J adr] + iudp = ap + ivp — W((J, v) + c, 

where the function W((T, v) is expressed in terms of the dilogarithm Li 2 {z), 
z/) = Li 2 {-e^^^G+v-q)-^ + ^ ^ 2 ^ 2 ^ 

Taking into account yet another classical formula. 


= —27ri\nG{z) — 27rz^ln 


sin(7r2;) 


TT 


TT 


-TT Z{1- Z) + y. 


where 


G{z) = 


G{l + z) 


G{l-zy 
and the elementary relation. 


( 68 ) 


(59) 


2cosir(<T + r(± j) = -‘f), 


we arrive at the following final expression for the first integral in (55) 

' G(a + q+i=p' 


(60) 


-Svri / adp + iudp = —Siiiap + 2 In 


G{a + p + ^) 


Anp — iTTU^ + 2i ln(27r)z/ 




27ia^ + QTrip^ + Aniap — 27Tia + 27Tip. 


(61) 


Substituting formulae (56), (57), and (61) in (55) we arrive at the equation. 


^ = ci(27rr2 


wq2i^^+(t^ 24—12(Tg27ri(?7^—2 (t?7—(T^+2?7—(t) 
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X 


r(l - 2a) ( G(1 + iu)G{l + 2a)G(l - 2a)G(a + r]+ 
r(2a) I G(a + r/ + 2^) 


(62) 


1 ^ 

where ci is a numerical constant. We know, that ifM = 0,a = r 7 = = 0, then r = const • 

and Goo = Gq. This choice of parameters satishes conditions (g. Hence, 


220- 


Cl = 


vr(G(i))4- 


To complete the proof of Theorem we only need now to lift the technical assumption (44). 
This can be justihed by noticing that the both sides of (62) are analytic functions of the 


Riemann-Hilbert data. (For the left hand side it follows from the general Birkhoff-Grothendieck- 
Malgrange theory.) 

Remark 5. The variables (p, a) and {—ip, v) are canonical variables. In fact, one has that 

U3. 

G = Z2nidr] A da = 327idp A du. 

The function W was introduced in JT^ as the generating function of the canonical transforma¬ 
tion 

(hW) ^ {-^P, ^)- 


Indeed, using (59), (60) and the fact that Li' 2 {z) = —z ^ \\i{z — 1), one can show that fWf, 


dW , . dW 

p = ——, and ip = 


da ’ 


dr 


The last eguation is also eguivalent to the integral formula 

In [12], and in fact earlier in the pioneering works [5], [6], the derivation of the constant 
terms in the asymptotics of the tan-functions was based on the heuristic assumption (followed 
from the conformal block representation of the tan-functions) that these constants are related to 
the generating functions of the relevant canonical transformations between the canonical pairs 
associated with different critical points. In the case of equation ([^ the points are 0 and oo and 
the generating function is the function W. This is a very important conceptual point, and our 
analysis justihes it in the case of the Painleve III equation ([^. It is also worth noticing that 
this hamiltonian interpretation of the ratio Goo/Go is already present in formula (54). Indeed, 
this formula tell us that the logarithm of the ratio Goo/Gq is, up to the elementary function, 
12 ^ -|- a‘^/8 — ir, the generating function of the canonical transformation between the Cauchy 
data {a, (I) and asymptotic at inhnity data (6+,5_). 


6 Proof of the ILT-Conjecture 

In [121 different r-function was introduced 


1 1 iu(x) 


Tra{2 ^V) = (r(a;))2x4e 4 , 


(63) 
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(1 + o(l)), a; —)■ 0, 


(64) 


r^(2-^V) = 


2 ” 


G{1 + 2a)G{l - 2a) 


Tra {2 = x(o', ^)e '4 2 ^"( 27 r) 2"G(1 + iz/)x'2+ 4 ei 6 +'^*(l + 0(1))^ x ^ 00. 


But from formulae (|ll[), 
7-^(2 


—12^4^ _ 2 ^—^+i7Tr}<2^ — Qa 


X ) = GqC 4 


we also have 
2r(l -2cr) 


X 


r(2a) 


(1 + 0(1)), X —)■ 0, 


(65) 


( 66 ) 


r™(2-'V) = C'.|x^+3ew+""(l + o(l)), x ^ cx). 


From formulae (64), (65), (66), (67) we get 
1 


C'.l ( 27 r)f 2 -i-"^-i 2 a 2 + 6 -e-^--^+x / r(2cr) 

^ G'(l + w)G'(l + 2cr)G(l-2a) Vr(l -2a) 


X(a, u,r]) = 


(67) 


Substituting here expression for we get the formula conjectured in [12] 

^ ’ G\\)G{a + r^+^^)^ 

7 Hamiltonian meaning of Malgrange-Bertola Differen¬ 
tial form 


One can rewrite (42) 


ui = 


xdx T-idx xT-Lpdp vupdp xT-Lqdq vugdq 


xdx d^xl-i) ^ xTixdx vdu ^ vUxdx 


It follows from ([^ that 


Using this formula we get 


4 4 4 

xT-Lx + vux = T-L. 


1 fx^ \ 1 

u = —d I-h xV, I H— ( T-Ldx — vdu 

4 V 2 / 4 


( 68 ) 


We want to emphasize that all the objects are considered as the functions of the triple (x,p, q) 
and all the differentials are taken with respect to all these three variables. 


From (|6^ it follows that up to the multiplication by -4 and the subtraction of a total 

(69) 


differential, the Malgrange-Bertola form u coincides with the form 
vdu — 1-Ldx = vUxdx + vUpdp + vuqdq — Tidx. 
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The restriction of this form on a trajectory of the Hamiltonian system Q, i.e. on the curve, 
p = const, q = const, (70) 

in the extended space of the monodromy data {{x,p, q)} coincides with the form 


dS{x) = vUxdx — Tidx, 


where S{x) is the classical action evaluated on the trajectory (70). Hence, the Malgrange- 
Bertola form uj can be treated as a natural extension of the canonical form T-Ldx — vu^dx. 
It follows then, that the tan-function itself can be identihed with the classical action. More 
precisely, along any classical trajectory, we have that 


dlnr = 


dlnr 
—;—dx 
dx 


1 

4 


dS 


-d(x'H) 


_ ^ d{xn) \ 

\ A dx A dx ) 


(71) 


Of course, this differential identity can be easily (after it is written) checked directly. In its 
turn, it allows us to write the following representation for the ratio C^jC^ in terms of the 
regularized action integral. 


In 


Oo 


lim lim 

^ 0—^0 —^~ l “00 


[I — 


io 


_L 

8 


— 2uti — Inti — 



(72) 


It is worth noticing that, unlike the integral J Hdx, the action integral suits well to the differ¬ 
entiation with respect to p and g; indeed, after the relevant integration by part the remaining 
integral term would disappear in view of ([^. Therefore, equation (72) provides us with the 
possibility of an alternative derivation of our key formula (53). This derivation would be very 


similar to the evaluation of the action integral of the McCoy-Tracy-Wu solution of the PHI 
equation in 


Remark 6. Observe that the extended (with respect to x, p, q) differential of the form vdu—Hdx 
is the symplectic form H. Therefore, the fact that the Malgrange-Bertola form differs from the 
form —\{ydu — Hdx) by a total differential is a fact of general theory; indeed, the (extended) 
differentials of the both forms coincide - they both are the same 2-form, i.e. — ^fl. The additional 
information we are obtaining in (68) is the explicit evaluation of this total differential. This 


allows us to relate the tau-function and the action differential explicitly which would be important 
for the alterantive evaluation of the tau-constant via the action integral. 


It might seem quite surprising that one needed to start with the Malgrange-Bertola form 
in order to discover a rather simple differential identity (71). The absence of the very idea 
that the logarithm of the tau-function might differ from the classical action just by a total 
differential partially explains this. We now believe that the similar fact should be true for any 
isomonodrtomy tau-function, although it has been apparently missing in the general monodrmy 
theory of linear systems 0 


^In 2000, the first co-author together with Percy Deift tried to use technique of [T31 for evolution of the 
constant factors in the asymptotics of the Painleve V tau-function associated with the sine-kernel. We failed 
then because we did not have the analog of the relation 0 for the Painlevve V tau - function we were working 
with. Perhaps, it would make sense to revisit the issue now (although, the relevant constant factors have already 
been evaluated since then). 
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8 Appendix. The Derivation of Estimates (46) and (47). 


Substitution of the extended asymptotics (45) into the right hand side of (42) leads to the 


following expressions for its individual terms. 

,2 U‘2 „2ix^2iu 


xul 


^2 ^—2ix^—2iu 


^2 ^—2ix^—2iu—l 

32 


y2 ^ix^2iv—l 


32 

,2 


+ 2z/ + 3i) 


(6w^—2i/+3i)—1/-) 


X 


64 


+ 


^4 ^—‘iix^—Aiv—1 


64 


fO(a: 


—2+6|Si/K. 


); 


X 


-(COSM — 1) = 
4 


^2 2ixx 


2iXr^2iu 


+ 


^2 ^—2ix^—2iu 




‘iix ^Aiv—1 


^4 ^—Aix^—Aiu—l 


lJ2 ^2ix ^2ii'—l 


X 


32 


(6w^ + 2v — i) 


^2 ^—2ix^—2iu—l 


32 


64 64 

(6^2 -2i/-i) + 


-UpSinu= -- 


9 ^ -N b_b+„x 

+ bv — i) ^ --dG 


+ 




16 


Sb^b+pC ‘^'■^x ^ 6_6+pZ/ ^ Inx ib'^Upe ^*^x ^*^+^lnx 

64 ^ 4 ^ 4 4 

+ ^ ^ -(6z/2+2iz/-l)—-(6z/2-2ii/-l) 


16 


16 


16 


b_pb+x ^ 6_6_pe-2“x-2*'^+i 36^6_pe2“x2*'^ 6^ 5_pe-^“x-^*^ ^ 6+6_pZ/ 


64 


16 


^2ix^2ii' 


ib^ UpC ‘^*'^lnx 6_6_pe ^*^x , o , b\vpe"'“"x" . . 

+ -7x--(6fz/2 - 5zy - i) + -(3iz/ - 1) 


16 


16 


b‘^Vpe "‘"x 


—2ix^—2iu 


.(3,Zy+l)+Zy^Zy + 0(x-l+6M). 


X 

T 


Uq sm u = ^ p —t g (>; 
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^—2ix^—2iu+l 


z6iz/pe^*^a;^*'^'''Mnx 6_6_p 

~^UxUpx 


+ 


4 4 

ib‘^i>pe~‘^^^x~‘^^^~^^ In a: 


p^ /o- l1^ b^b^pC^^^x'^'"' 2 Q-\ b_b+pX 


(Siu+l)- 


-{Qii'‘‘+iy+3i) + 


16 .' 4 

Inx, , , b^b+pC 

+ + ^ _ -(6i/2-2iz/ + 3) + ^ 


bj^b_ph> b^b+pV 6^i/pe^“x^*^lnx^^ 2 


b^b_pX 

- ~\ -1- 

4^44 16 

Inx 6^ Inx, , x &-&-pe“^“x“^*^ , n 

■ ^ -(6z/2+2iz/+3) +-^ 


42X^4^^' 


16 

-, — 22X^ — 221/ 


16 


16 


16 

{6iu^ — u+3i) 


6^ UpC ^®*x , , b^ 6_pe ’""x 

~ ^ ^ -(3zz/ - 1) + 


i6iz/pe-4“x-4*^lnx 


16 


16 


b^b^pC 


— 2ix^—2ii/ 
Jb 


2ix^2ii/ 


64 


b^b-pC^-'x 

U 


+ 0(x-^+®l®"l); 


X 

4 


— UxUqx = \ p^ q\] 


xUxUp ifo+fo+pC^^'^x^*^ 6iz/pe^**x^*^ Inx ib+b^p , i6_6+p 

—-^-^-“ 2z/pZ/ In x-^ 

4 4 4 4 ^4 

i6_6_pe“^*®x“^*^ 6^z/pe“^“x“^*^lnx ^ _i+4|Qp|^ 

;; ^ h 0(x ); 


XUxUq 


= iP^Qh 


It is quite remarkable that the substitution of these long expressions into the right hand side 


of (42) yields to the compact formula (46). 


For asymptotics at zero we get the following estimates. 


a 


g ^(cosu-1) = ^ + 0(x^ 


X^UpSinU ^ 3 ;% Sinu ^ ^^^2-|S(o 




X^UxUpx _ Otap ^ 2-|Q(«)K x'^UxUqx _ OtOiq , 2-|S(«) 


X 


'i 

xUxUp aap In X ^ a/3j 


X 




4 4 4 

^'^xUq ^ aaq InX ^ ajbq ^ 2-|S>(a)| 

4 4 4 ^ ' 


These equations yield at once (47) 
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